Fisher information and spin squeezing in the Lipkin-Meshkov-Glick Model 
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Fisher information, lies at the heart of parameter estimation theory, was recently found to have 
a close relation with multipartite entanglement (Pezze and Smerzi, Phys. Rev. Lett. 102, 100401). 
We use Fisher information to distinguish and characterize behaviors of ground state of the Lipkin- 
Meskhov-Glick model, which displays a second-order quantum phase transition between the broken 
and symmetric phases. We find that the parameter sensitivity of the system attains the Heisenberg 
limit in the broken phase, while it is just around the shot-noise limit in the symmetric phase. 
Based on parameter estimation, Fisher information provides us an approach to the quantum phase 
transition. 



I. INTRODUCTION 

Parameter estimation of probability distributions is 
one of the most basic tasks in information theory, and 
has been generalized to quantum regime 

[El 

since the 

description of quantum mechanics is essentially proba- 
bilistic. How to improve the precision of parameter es- 
timation has been focused for many years, and is of im- 
portant applications in quantum technology like quan- 
tum frequency standards 0, E| , measurement of gravity 
accelerations Q, and clock synchronization [|| etc. 

Consider a quantum state pg = U$p- lu Ug, where Ug = 

cxp (iOK^j and K is a generator. We estimate parameter 

through proper measurements, however, the precision 
of our estimation is limited by the quantum Cramer- Rao 
(QCR) bound 



In the past, many works have been devoted to improve- 
ment of parameter sensitivit y by using entangled states 

IiiIiilSliSISEIiil^l2il^l2al2j. 

Quite recently, Pezze and Smerzi [25| found an inter- 
esting application of QFI in multipartite entanglement 
and the sub-shot-noise phase sensitivity in the estima- 
tion of a collective rotation angle. Consider an ensemble 
of spin-half particles in the state p ln , they introduced a 
quantity 



N 



F(pin,Sft) 



(2) 



A# > (A0) QCR = 



vF ( /?i n , K 



(1) 



and prove that % 2 < 1 implies multipartite entangle- 
ment. Here, the generator of 9 is = S ■ n that 
denotes the collective spin operator along direction n. 
Namely, a sufficient condition is given for quantum en- 
tanglement. We may define a mean Fisher information 
as F m = F (pi n , Sft) /N. Then, \ 2 an d F m are reciprocal 
to each other. The relation between \ an d QCR bound 
is 



9 is the so-called unbi- 
is the quantum Fisher 



where v is the number of trails 

ased estimator, and F{p m ,K\ 

information (QFI) [l], @, 0, S|- In a sense, parame- 
ter estimation is equivalent to distinguishing neighbor- 
ing states along the path in parameter space. We know 
QFI has close relation with Bures distance @ , the most 
studied distance in quantum-state space, and Bures dis- 
tance is directly related to the Uhlmann fidelity [lfJl. 
The QFI is proportional to the Bures distance [IJ, LL2] • 
For pure states, the QFI, as well as the Bures distance 
ds 2 B , is just proportional to the variance of K 8], that 

is F(p m ,K) = Ads% = A^Aky 
increasing experimental times v, we can improve the es- 
timation precision A9 by choosing proper states p m for 
a given K. In general, entangled states are more sensi- 
tive than separable states, i.e., the variance of K is large. 



(3) 



N is the shot-noise limit and we set 

2 



where (A#) gN 

v = 1. Thus, it is evident that \ z < 1 becomes a nec 
essary and sufficient condition for sub-shot-noise phase 
estimation. 

In this work, we study the Fisher information of the 
ground state of the Lipkin-Meshkov-Glick (LMG) model 
26] , which has a second-order quantum phase transition 
(QPT) (27[, between a symmetric (polarized, h > 1) 
phase and a broken (collective, h < 1) phase. Some 
works have been devoted to study the LMG model with 

mm 



Therefore, besides concurrence [28( and entropy 
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In our work we 
find that, besides indicating the critical point and entan- 
glement, x 2 reflects the performances of ground states 
of these two phases in the sense of parameter sensitiv- 
ity. In the symmetric phase, \ 2 approaches to 1 with the 
increasing of h, and is independent of N, that means 
(A#)q CR ~ (A0) gN . In the broken phase, we find 
X 2 — thus (A#)q CR ~ 1/N attaining the Heisen- 

berg limit. 
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This paper is organized as follows. In Sec. II, we 
give brief discussions about the relations between spin 
squeezing and \ 2 . Then in Sec III, we study \ 2 and 
spin squeezing for the ground state of the LMG model in 
both isotropic (7 = 1) and anisotropic (7 7^ 1) cases. In 
isotropic case, the LMG model is diagonal in Dicke states. 
For Dicke states, we find that, % 2 an d the spin squeez- 
ing parameter by Kitagawa and Ueda [3l| are recipro- 
cal to each other. In anisotropic case, we use Holstein- 
Primakoff transformation and derive % 2 in the thermo- 
dynamic limit. The finite size behaviors of y 2 and the 
spin squeezing parameter in the critical point are also 
obtained. The numerical results coincide well with the 
analytical ones. 



II. FISHER INFORMATION AND SPIN 
SQUEEZING PARAMETERS 

Fisher information is related to spin squeezing, and 
there are two spin squeezin g p arameters respectively 
given by Kitagawa and Ueda |31| , and Wineland [12] , 



2 _ 4(A5^ 
4l _ N 



N(AS H± f 

\(SnW 



(4) 



where subscript n± refers to an arbitrary axis perpendic- 
ular to the mean spin (S), where the minimum value of 
(AS*) 2 is obtained. The inequality £ 2 < 1 (i = 1,2) in- 
dicates that the state is spin squeezed. Spin squeezed 
states can be used to reduce the measurement uncer- 
tainty [3lL [33| , and improve the measurement precision 
of the atomic clock transition [34], |35| . The spin squeez- 
ing inequality is a criteria for multipartite entanglement 
[39. [37| . For an arbitrary multiqubit separable states, it 
was found that £f > 1, and thus £f < 1 implies quantum 
entanglement. 
As proved in (25j . 

F(p,5 sl )(A5 s J 2 >|(5 s )| 2 , (5) 

where directions n'^, ri±, n are orthogonal to each other. 
F ^p, S.fi> ± ^ = 4 ( ARj , where R is determined by Rp + 

pR = 1 (j>Sit' ± ~ S ff± p). In general, (ar) < (ASr.J 2 , 
and the equality is obtained only for pure states. Then 
Eq. ([5|) reduces to the usual uncertainty relation, 



(AS H1 ) 2 (ASn ± Y > 



\(Sn) 



(6) 



for pure states. The above inequality can be written in 
terms of the inverse of the mean QFI and the squeezing 
parameter £f as 



N(AS. n± Y 



> 



N 



4(A%, ) 2 



^ = X 2 - IT) 



Both the inequalities, £| < 1 and mean QFI F m > 1 
(x 2 < 1), indicate the presence of entanglement. 

Furthermore, £ 2 < £ 2 , and there is no similar relation 
between £ 2 and x 2 like Eq. (|7|). However, we find that 



£x 2 



(AS Hj 



< 1, 



(8) 



since (AS*^^) 2 ((ASjt ± ) 2 ) is the maximum (minimum) 
variance. As proved in [37| . if the pure state is of ex- 
change symmetry, £ 2 < 1 implies entanglement. Then, 
from the above equation % 2 > 1 implies entanglement. 
We know that \ 2 < 1 indicates entanglement too. There- 
fore, a pure symmetric state is entangled iff % 2 7^ 1 

When the mean spin direction is along z direction, the 
squeezing parameter £ 2 and \ 2 becomes 



Si N 



N 



4 max (S 2 



(9) 



where 



S± = coseS x +sm6S v . (10) 



Furthermore, if ({S x , S y }) = 0, for instance, in the LMG 
model [28j, we have 



, 2 _ 4min((Sj),(Sg)) 2 
?i N , X 



N 



4max((5 2 ),(^ 2 )) 



(11) 

thus we only need to compute (S 2 ) and (S 2 ) to deter- 
mined the squeezing parameter and quantity \ 2 i n t ne 
following discussions of QPTs in LMG model. 



III. FISHER INFORMATION AND SQUEEZING 
IN THE LMG MODEL 



The LMG model, originally introduced in nuclear 
physics and has found applications in a broad range of 
other top ics: statistical mechanics of quantum spin sys- 
tem [23], Bose-Einstein condensates [39|], or magnetic 
molecules such as Mni 2 acetate Recently, some 

quantum-information concepts, such as q uantum entan- 
glement [29|, [lOl and quantum fidelity |4ll.l4^|. have been 
studied in this model, aiming at characterizing its QPT. 
It is an exactly solvable [43|, |44[ many-body interacting 
quantum system as well as one of the simplest to show 
a quantum transition in the regime of strong coupling. 
The quantum phase transition of this model is also clear 
and interest: the ground state becomes degenerate and 
a macroscopic change in the ground state energy takes 
place. 
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LMG Hamiltonian 



The Hamiltonian of the LMG model reads 



H 



1 

N 



(S 2 + 7S 2 ) - hS z , 



(12) 



where S a = X^iLi a a/^ are the total spin operators in 
the direction a — x,y, z; a % a are the Pauli matrices; N is 
the total spin number, 7 is the anisotropic parameter and 
h is the effective strength of the external field. Without 
loss of generality, we assume < 7 < 1 and h > 0. 

This system undergoes a second-order QPT at h = 1, 
between a symmetric (h > 1) phase and a broken (h < 1) 
phase, which are associated with single-particle and col- 
lective behaviors, respectively. These two phases are well 
described by a mean- field approach [28]. The classical 
state is fully polarized in the field direction ((cr*) = l) 
for h > 1, and is twofold degenerate with (a),} = h for 
h < 1. However, since the Hamiltonian is of spin- flip 



symmetry, i.e., 



A' 



= 0, we have 



(S X ) = (Sy) = 0, (S X S Z ) = (SyS Z ) = 0. 



(13) 



Thus the mean spin direction is along the z-axis for finite 
size case. In addition, \H, S 2 ] = 0, and the ground state 
lies in the S = N/2 symmetric section. 



B. Isotropic case and Dicke state 

We begin with the simple isotropic case, 7=1. The 
Hamiltonian reduces to 



H 



1 

'N 



(S 2 - S 2 ) - hS z 



(14) 



which is diagonal in the standard eigenbasis {\S, M}} of 
S 2 and S z . For S = N/2 the energy eigenvalue is 



(«-f) S -f( 1 + " 2 )' < 15 > 



and the ground state \S,Mq) is readily obtained when 

m 

f N/2 for h>l, 

M ° = \N/2- R[N(l-h)/2] for < h< 1, ^ 16 ^ 

where R(x) = round(x) gives the nearest integer of x. 
Then one can see level crossings exist at h = hj, where 
hj = 1 — (2j + 1) /N, between the two states \S, S — j) 
and \S,S-j-l). 

As the ground state is actually a Dicke state \S,M), 
(S 2 ) = (S 2 ) = (S 2 + S-M 2 )/2,then 

N 



2(S 2 + S~M 2 ) 
1 



N/2 + 1- M 2 /N 



< 1, 



(17) 



the equality is obtained for M 
have 



±S. Immediately, we 



e? - 1/x 2 > 1. 

As we know that, when M 7^ ±S, the Dicke states are 
entangled but not spin squeezed, since £f > £ 2 > 1. Nu- 
merical results of £ 2 and \ 2 f° r the isotropic LMG model 
are shown in Fig. [T] (d). We can see that, in the broken 
phase, M < S, £ 2 = 1/x 2 , while in the symmetric phase, 
the ground state is \S, S), thus x 2 — £1 = 1- 

By considering % 2 in Eq. (jTTj) . when M is close to ±5, 
X 2 is just a bit lower than 1, thus A9 is not improved 
much than (A0) SN . When \M\ < y/S, we have 



< x 2 < 2/ (JV + 1) , 



(18) 



and thus 



(A0) 



QCR 



which attains the Heisenberg limit. Although |S, ±S) is 
not entangled, "cat state" (or GHZ state) 



= (\S,S) + \S,-S))/V2 



(19) 



is entangled and useful in phase estimation [161 ] . Un- 
der (S a ) = 0, for a = x,y,z, thus there is no spin 
squeezing. We find the maximum variance (AS Z ) 2 = S 2 , 
then 



X 2 W = 1/N, (A8) QCR = 1/N, 



(20) 



beating the Heisenberg limit. From the above analysis 
we know that, for typical symmetry multipartite states, 
Dicke states, there are no spin squeezing, while % 2 < 1 
indicates that they are entangled and are useful resources 
for phase estimation. 



C. Anisotropic case 

Now we consider the anisotropic case, < 7 < 1. The 
spin expectation values (S 2 /) can not be obtained ana- 
lytically. By treating the quantum effect as small fluc- 
tuations, approximate results can be obtaine d by using 
the Holstein-Primakoff (H-P) transformation [45|] in the 
thermodynamic limit, and b y u s ing the continues unitary 
transformation method [4611471 l48j for finite size case. 

In the thermodynamic limit, the quantum fluctuations 
are small, we can use the H-P approximation. This 
method requires one to determine the semiclassical mag- 
netization (S), which is not along z-axis in the broken 
phase in the thermodynamic limit. Following conven- 
tional steps, we first employ a mean field approach, define 
a spin coherent state 



N 

= <g> 



cos -|0) ( + e 



i<t>/2 



sin-|l)j 



(21) 
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under which 



N 



(9,4>\S\6,<p) = — (sin^ cos <f>, sin0 sin <f>, cos 6). (22) 



The Hamiltonian is rewritten as 



H = -N 



^ sin 2 6 (cos 2 ( 



7 sin 



h cos 6 



(23) 



As max [cos 2 <f> + 7 sin 2 ^] = max (1, 7) = 1, (7 < 1), we 
have 



mini? = — TV max 



■ sin 8 + h cos 9 



(24) 



then we conclude: (i), symmetric phase, h > 1, #0 = 0, 
for all 7; (ii), broken phase, < h < 1, #o = arccos/i, 
= 0, 7T, for 7^1. We emphasize that, the mean spin 
direction is along the z-axis when the system is finite. 

We rotate the z-axis to the semiclassical magnetiza- 
tion, 



Sy 

1 Sz 



cos $0 sin#o\ / S x ^ 



Sy I = 



1 



S 1, 



sin6» cos6» / \§, 



(25) 





10' 



10" 



10" 



(d), y =T% 


— * 2 











FIG. 1: and x 2 as functions of h for various 7, with system 
size N = 100. The crossing points of £1 and the horizontal 
line in the broken phase are h — ^fy. 



As presented in [28], 0q = for ft, > 1 so that S = S, and 
6q = arccos/i for h < 1. The transformed Hamiltonian 
reads 



H = — hmS z — 



1 



/iVT 



™ 2 +7 52 3m 2 + 7 - 2 - 
S 2 

m 2 — 7 



5+ + 5_ - 



AN 



si + s*_ 



1 vi 



27V 



S+S z + S Z S+ + S~S Z + S Z S~ 



(26) 



where m = cos#o- Then we introduce the H-P transfor- 
mation 



S z =N/2-a ! a, S- = VNaty/l - ata/N = 



t 

(27) 



The Hamiltonian can be written as 



£ (0) = 2ftm-3m -7 + 2^ _ m_7 (flta fl2) 
2 4 v ' 



+ 



1 — m 2 



(28) 



up to the 0-th order of AT. We neglect the terms of the 
1-th order of N as they are constant. Now we use the 
Bogoliubov transformation 



a) = cosh ( - ) 6* + sinh ( - ) b. 

2 7 V2/ 



To diagonalize £T (°) , we find 



tanh t 



m — 7 



2hm — 3m 2 — 7 + 2 



(29) 



(30) 



The rotated spins are written under the H-P represen- 
tation, 



y/N (\ +e 



Si, — 



\ 1/4 

2 \l-e) 



2 \l + e 



(6t + b)+0 (l/N) , 
(ft* + 6) + O (1/iV) , 



TV 1 ^ 1 



22V VT 



1 



6 t 6 + £ ( 6 t2 + ^ 



For symmetric phase, m = 1, S a = S a , we have 



<s2> 
<^ 2 > 



5't 



S? 





7 


h.- 


1 


/l- 


1 


/l- 


7 



(31) 



(32) 



while for broken phase, m = h, we need to rotate back 
to S x as 



^ = cos 9q + S z sin 0q 



h~S x + - h?S x 



(33) 
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FIG. 2: Comparison between analytical and numerical results 
for £1 and \ 2 , in the case of 7 = 1/2. The solid line is ana- 
lytical results for the thermodynamic limit (TL). 



then wc have 

(s 2 x )=VT^(s 2 ) + h(s 2 x ) 

TV (1 - 7) ft 2 - (2 - ft 2 - 7) (1 - h 2 ) 



^(l-ft 2 )(l- 7 ) 



(34) 



We insert the above results into Eq. (fTTj) . For polarized 
phase, ft > 1, we have 



t-2 2 

£1 = X 



ft — 7 



< 1. 



(35) 



When ft, is far from the critical point, £f and x 2 approach 
to 1, then A# ~ (A6*) gN . For broken phase, ft < 1, we 
get the spin squeezing, 



while 



N 



1-ft 2 
1-7 ' 



1 



A{Sl) ~ {N + 2){l-h?) 



(36) 



at ^ 



Thus (A0)q CR ~ 1/N. When ft approaches to the crit- 
ical point ft c = 1, there are two limit processes in Eq. 



10" 



10 



10 



10" 
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N=£(JU 
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---N=300 
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0.5 



1.5 



FIG. 3: Plot of x 2 as a function of h with 7 = 1/2 for different 
system sizes N = 100, 200, 300, 400, (from top to bottom). \ 2 
is nearly independent of JV for h > 1, while it drops with the 
increasing of N in the broken phase, h < 1. 



that is (1 — ft) tends to be zero and TV tends to be 
infinity. To overcome this problem, we need to expand 
the Hamiltonian in higher order of 1/N. Fortunately, the 
finite size behaviors of the spin squeezing and x 2 a t the 
critical point can be derived by using the results obtained 
in [28], where the authors employ the continues unitary 
transformations and get 



4(S 2 



N 2 



a {0) 



h=l 



TV 2 / 3 ' N 2 



h=l 



a {0) 
a yy 



(38) 



where a^/1 and dyy are constant independent of N . Now 
we have 



e 2 \ 



a (0) 
u yy 

N 2 / 3 ' 



1/1=1 



a (0) 

AT2/3' 



(39) 



then for large AT, £ 2 and x 2 converge to zero as 1/iV 2 / 3 , 
and (A0) QCR ~ 1/7V 5 / 6 . 

To verify these analytical prediction, in Figs. Q] and [21 
we show numerical results for £ 2 and x 2 as functions of ft 
with different 7 for finite size system. As shown in Fig. [TJ 
in the symmetric phase, £ 2 = x 2 < lj while in the broken 
phase, x 2 an d behaves very differently. In the most of 
parameter range, % 2 < 1, which indicates entanglement, 
while for ft < ^7, £ 2 > 1, and thus the system is not spin 
squeezed (Fig. 1(d)). For the isotropic case, there is no 
spin squeezing. In Fig. fJl we plot x 2 an d £ 2 for N = 500 
and the thermaldynamical limit, and find the numerical 
results coincide well with the analytical ones obtained by 
H-P transformation method. 

As shown in Figs. [3] and [3J x 2 is nearly independent 
of larger N in the symmetric phase, and approaches to 
1/N as ft being away from the critical point in the bro- 
ken phase. Therefore, entanglement characterized by x 2 
is very different in the two phases, especially when we 
treat them as resources for quantum estimation. In the 
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FIG. 4: Scaling of 1/x 2 as a function of N with h = 7 = 1/2. 
The diamonds are numerical results, while the solid line is a 
linear fit with a a constant. It's obviously that x 2 c< 1/iV. 



symmetric phase, as shown in Fig. [3J \ 2 i s nearly inde- 
pendent of N, and the parameter sensitivity is at the level 
of (A0) gN , while in the broken phase, the ground states 
are more sensitive in parameter. In Fig. ([4]), we show 
numerical results for x 2 m the broken phase at h = 1/2, 
7 = 1/2, and we see clearly that x 2 oc 1/N. Therefore, 
the parameter estimation in the broken phase is enhanced 
to the Heisenberg limit. 

One can use the concurrence and entropy to quantity 
entan glem ent in the LMG model and results are obtained 
in [28l I29tl3fj| . We see that, both concurrence and entropy 
indicate well the presence of entanglement, however, from 
them, we cannot tell whether the entanglement of the 
ground state is useful in parameter estimation. From re- 
sults of x 2 ) we can see that, the entanglements in these 
two phases are different according to their performances 
in estimation. On one hand, we can use the collective 
behavior of the LMG model to improve the phase esti- 



mation precision, on the other hand, the differences of 
the parameter sensitivities can be used to characterize 
and distinguish the two quantum phases. 



IV. CONCLUSION 

We have analyzed x 2 an d spin squeezing parameters £ 2 
in the ground state of the LMG model. For isotropic case, 
the Hamiltonian is diagonal in Dicke states, for which we 
have £ 2 > £ 2 = 1/x 2 - For anisotropic case, our results 
indicate that, x 2 classifies states in different phases in 
the sense of quantum phase estimation. Hence we can 
use x 2 to distinguish and characterize the behaviors of 
the two phases of the LMG model. In the symmetric 
phase, x 2 is independent of N and approaches to 1 with 
the increasing of h, thus (A0)q CR ~ (A6*) gN , that is just 
a bit lower than the shot-noise limit. In the broken phase, 
we find x 2 — V-W anc ^ (^)qcr — V-^j which attains 
the Heisenberg limit. 

Fisher information, being related to the Cramer-Rao 
inequality, is used to measure how much information that 
we know about some certain parameters in a probability 
distribution. From present results, we see that Fisher in- 
formation can also characterize the QPT, by distinguish- 
ing the entangled ground states in the sense of parameter 
sensitivity. This approach is promising and expected to 
be applicable to other spin systems undergoing a QPT. 
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